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ABSTRACT: The general assumption is that stiffness matrix is constant in structures. But if the loads 
increase, deformation-load curve becomes nonlinear. The aim of this study is to analyze the vibration of a 
structure under these facts. The chosen sample problem is a steel single storey frame with a horizontal force. 
At first, classical stiffness matrix is determined under static loads. Then the structure has been reduced to a 
unique mass-spring system with increasing loads. It is accepted that material is ideal elasto-plastic. The change 
of stiffness under increasing loads will be calculated by plastic analysis. The system is not like the beginning 
after first plastic hinge occurs and degree of freedom of the system changes. If the load increases the decreas- 
ing of stiffness matrix continues. It is accepted that stiffness matrix is constant until first plastic hinge occurs. 
Also, it is accepted that the stiffness matrix is also constant between first and the second plastic hinges. By this 
way, a curve is obtained between load and stiffness matrix and the stiffness-force diagram has been plotted. 
In addition, the positions of plastic hinges have been controlled by SAP2000 computer program. Finally, the 
forced vibration of the sample structure under an harmonic load has been investigated. It is clear that the 
external load varies by time. Then the force-time diagram and the stiffness-time graphic have also been plotted 
under increasing loads. The damping coefficient of the system must be calculated for every stiffness value by 
choosing a initial value for system damping. The displacement-time curve has also been given. The operations 
of increasing loads have been repeated for decreasing loads. Also the displacement-time curves for increasing 
and decreasing loads have been combined and they have been shown on a single diagram. 


1 INTRODUCTION tion is made using certain models. The main fac- 
tors that determine results, the distribution of 
The aim of this study is to calculate the variable the masses in the structure, the characteristics of 
stiffness matrix under increasing and decreasing external forces and displacements, internal fric- 
loads. As the first step the solution of any struc- _ tion, cracking, the resistances against deformation 
ture under static loads is solved by classical matrix- of the bars(elements) and damping mechanism. 
displacement method assuming stiffness matrix is _ (Yerlici and Lus, 2007). 
constant. And the stress resultants and stresses are Consider a system which can only rotate about 
calculated in the bars forming the structure. Then an axis or can move in just one direction. This 
the loads have been increased. As a result of this model is called a single degree of freedom system. 
increment a specific cross section of a certain bar The dynamic behavior of a structure under exter- 
cannot carry more moment. Then a plastic hinge _ nal influences that varies depending on the mass of 
occur at this point. The new system is different system, the stiffness and energy loss in the system. 
from the original system. Again load is increased To solve a system under dynamic loads the equa- 
and new system is solved till a second plastic hinge __tion of motion for the system is necessary. One 
occurs. Here only plane systems are considered and equation is sufficient to determine the dynamic 
the moment at a plastic hinge is bending moment. _ behavior of a single degree of freedom system. 
Degree of freedom of the first system decreases one 
degree after each plastic hinge occurs. Then a specific 
point is selected and the deformation of this pointis 2 SAMPLE PROBLEM 
calculated for every stiffness matrix. These procedure 
will continue till the system is unstable. And the load- — The chosen system is a plane frame with three bars 
deformation curve is plotted for this specific pointand — which is shown in the Figure 1. A, and A, ends 
system is reduced a mass-spring system. And using are connected to outer medium by built in connec- 
deformation-load curve the variation of the stiffness tions. The chosen material is steel with young mod- 
of this imaginary spring by load is obtained. ulus 2.1 x 10’ N/cm? and density 7850 kg/m? and 
After these same structure is examined under yield stress is 275 N/mm”. The material is assumed 
dynamic loads which vary by time. Dynamic solu- as ideal elasto-plastic. Degree of freedom is three. 
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Figure 1. 


Sample problem. 


Cross sections are square 6 x 6 cm for beam and 
8 x 8 cm for columns. 

There is only a horizontal P, force acting at 
point B,. Horizontal displacement U,, at point B, 
is determined by solving this system via matrix- 
displacement method. Here only this displacement 
is considered then the stiffness of the system is 
defined as 


(1) 


Then the P load is increased and at one P, value 
first plastic hinge occurs. It is assumed that stiff- 
ness is constant till this value of P, First hinge is 
at A, point. 

This system is also solved and new stiffness k, 
becomes 


(2) 


Then this procedure continues till last plastic 
two hinges occur at A, and B, At this time this 
value of P is named as limit load. The variation of 
k versus P is given in Figure 2. 

Then the dynamic loading is considered for the 
spring-mass system given by Figure 3 

For this system equation of motion is 
mx +cx+kx = P(t) (3) 

Here k is the stiffness, m is the mass, c is the 
damping coefficient of the system and the force is 
assumed as a harmonic load as follows 


P(t) = Pcos(Qt —-@) = Poos(==1 —@) 


F 


(4) 


Here Q is the frequency, T;,is period and © is 
the phase angle which does not make any effect on 
reactions of the system. Then this term is ignored. 
The natural vibration of the system is also ignored. 
The chosen T,, is 0.4 sec for this problem. 
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Figure 2. The variation of stiffness versus force. 
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Figure 3. 


Spring-mass system. 


Two new variables defined 


(5) 


(6) 


Here € is damping ratio and @ is natural 
frequency. 

Three kinds of solutions can be found for equa- 
tion (3) these are 


e €>1 heavy damping case 
e €=1 critical damping case 
e & <1 weak damping case 
Here weak damping is considered and the solu- 
tion of equation (3) becomes 


a, =0(-gFiVi-#) (7) 
a=-€0 
b=ai-2 (8) 
x=e" (Guoe cosbtdt+C, jin btt+ 

(9) 


(-F1 P(the“sinbtdt+C, Jos bt 
0 


It is clear that the external load varies by time 
in the dynamic loading. Therefore, the time val- 
ues have been evaluated for increasing load values 
using the load values at which the plastic hinges 
occur. 
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Mass of the system is taken as 


m=115.23 N.s?/m (10) 

Then the force-time diagram has been plotted 
under increasing loads in Figure 4. 

After that, the stiffness-time graphic has also 
been given which shows time intervals of this 
stiffness. 

Natural frequency values of the system have 
been calculated for every time interval. 


Ke Q; > @, = 49.698 rad | sec 


(11) 
m 
= = @ > w, = 49.668 rad / sec (12) 
m 
ae 
5 = @ > w, = 49.659 rad | sec (13) 
m 


The damping coefficient of the system must be 
calculated for every stiffness value by choosing a 
initial ¢ value for system damping. 


& = 0.02 (14) 
Table 1. The variation of load versus time. 
Load (N) Time(sec) 
P, 33210 0.095 
P, 33225 0.096 
Ps 33300 0.1 
PSinét(N) 
5000 
30000 
25000 
20000 
15000 
10000 
5000 
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Figure 4. The variation of force versus time. 
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Figure 5. The variation of stiffness versus time. 
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Figure 6. The variation of natural frequency versus 
time. 
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Figure 7. The variation of damping ratio versus time. 
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Figure 8. The variation of displacement versus time. 
~ = 2£@, > ¢=229.068N.s/m (15) 
m 
~ =2E@, > € =0.0200121 (16) 
m 
~ ~2€, > & =0.0200157 (17) 


m 


According to this, the natural frequency-time 
and damping ratio of system-time graphics have 
been plotted. 

It is assumed that the stiffness matrix remains 
constant between two values of the load when 
using the stiffness-force diagram. The system has 
been solved in a time interval that corresponds to 
these loads. The values, at the end of the previ- 
ous interval have been used as initial conditions to 
determine the integration constants of the follow- 
ing interval. And the displacement-time curve has 
been plotted in Figure 8. 


Insights and Innovations in Structural Engineering, Mechanics and Computation — Zingoni (Ed.) 
© 2016 Taylor & Francis Group, London, ISBN 978-1-138-02927-9 


Table 2. The variation of load versus time. 


Load(N) Time(sec) 
P, 33300 0.104 
P, 33225 0.105 
P, 33210 0.2 
PSin& tN) 
5000 
30000 
25000 
20000 
15000 
10000 
5000 
0.05 0.10 015 on “°° 
Figure 9. The variation of load versus time. 
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Figure 10. The variation of displacement versus time. 


The operations of increasing loads have been 
repeated for decreasing loads. The changes of 
increased and decreased loads depending on time, 
has been represented in a single curve. 

During calculations, the stiffness values of every 
time range that is calculated for increasing loads 
and decreasing loads have been taken as the same. 
The displacement-time curve, for decreasing loads 
has been plotted under these conditions. The dis- 
placement-time curves for increasing loads and 
decreasing loads have been combined and they 
have been shown on one diagram. 


3 RESULTS 


The variation of the stiffness matrix by increas- 
ing loads has been calculated. And this variation 
is used for dynamic loading. The selected problem 
has been reduced a single mass-spring system. But 
if a set of system deformations is considered, stiff- 
ness transforms to a matrix. Same procedure can 
be extended to all calculations. But a restriction is 
necessary. Whole exterior loads are horizontal and 
have the same time dependencies. 
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